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I. Introduction 

In 1949, Levinson established a theorem in nonrelativistic quantum mechanics[l]. The 
theorem gives a relation between bound states and scattering states in a given angular 
momentum channel /, i.e., the total number of bound states n; is related to the phase shift 
Si{k) at threshold (A; = 0): 

Si{0)=ni7r, Z = 0,1,2,.... 
The case / = should be modified as 

So{0) = (no + l/2)7r 

when there exists a zero-energy resonance (a half bound state) [2]. This is one of the 
most interesting and beautiful results in nonrelativistic quantum theory. The subject 
has been studied by many authors (some are listed in the Refs.[2-8]) and generalized to 
relativistic quantum mechanics[6,9-14]. However, most of these authors deal with the 
problem in ordinary three-dimensional space. A two-dimensional version of Levinson's 
theorem does not appear to have been discussed by previous authors. In view of the wide 
interest in lower-dimensional field theories in recent years, e.g., Chern-Simons theory in 
2-1-1 dimensions, and the previous applications of Levinson's theorem to field theories[15], 
it seems of interest to study the theorem in two-dimensional space. On the other hand, 
as the problem exhibits some new features in two spatial dimensions, it may also be of 
interest in its own right. We are thus led to consider the problem. 

In a recent work [16] we have established the Levinson theorem in two spatial dimen- 
sions, which takes the following form: 

^m(0)=n^7r, m = 0,1,2,... 
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where 7^^1(0) is the phase shift of the mth partial wave at threshold, and rim is the total 
number of bound states with angular momentum m (it also equals the total number of 
bound states with angular momentum —m when m ^ 0). As in three dimensions, the 
modulo-TT ambiguity in the definition of r]m{k) has been resolved by setting 77^^(00) = 
(rather than a multiple of tt) which can be freely done in nonrelativistic theory. The 
theorem is similar to the three-dimensional one but simpler in that the existence of half 
bound states (possible for m = 0, 1) does not alter the form of the theorem. 

In this paper we extend the previous work to the relativistic case and establish the 
Levinson theorem for Dirac particles moving in an external central field in two spatial 
dimensions. In a given angular momentum channel j (j= ±1/2, ±3/2,...), the theorem 
relates the total number of bound states rij to the phase shifts rij{±Ek) at zero momentum: 

Vj{l^)+Vj{-I^) =11^71. (1) 

In the relativistic theory, one is not allowed to set 77^ (±00) = 0. But the modulo-7r 
ambiguity in the definition of rjj{±Ek) may be appropriately resolved (see Sec. V). As 
in the nonrelativistic case, the theorem is similar to the three-dimensional one[10] but 
somewhat simpler. In three dimensions, the theorem should be modified when there 
exists a half bound state, but here we have no such trouble. 

Throughout this paper natural units where h — c — 1 are employed. In the next 
section we first discuss various aspects of the solutions of the Dirac equation in an external 
central field in two spatial dimensions. Then we give a brief formulation of the partial- 
wave method for potential scattering of Dirac particles. In Sec. Ill the behavior of the 
phase shifts near A; = is analyzed. In Sec. IV we estabhsh the Levinson theorem using 
the Green function method[3,4,6,10]. Sec. V is devoted to some discussions relevant to 
the theorem. 
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II. Dirac particles in an external central field in two dimensions 

We work in (2+l)-dimensional space-time. The Dirac equation in an external vector field 
Ai,{t, r) reads 

(z7^D,-/.)* = 0, (2) 

where /j, is the mass of the particle, = + ieA^, e is the coupling constant, and 
summation over the repeated Greek index i/ (i/ = 0, 1, 2) is imphed. The 7^ are Dirac 
matrices satisfying the Clifford algebra: 

{7^ 7^^} = 2<7^^ (3) 

where g^^ = diag(l, —1, —1) is the Minkowskian metric, fn this paper we only consider the 
zero component of A^, which is cylindrically symmetric, namely, we consider the special 
case where 

A = 0, eAo = y(r), (4) 

where r is one of the polar coordinates (r, 9) in two-dimensional space. In this case we 
may set 

*(t,r) = e-^V(r), (5) 
and get a stationary equation for ip{r): 

Hijj = Ei/j, (6) 

where the Hamiltonian 

// = Q;-p + 7V + 1^(r), (7) 

where p = — iV, a — 7^7 or a* — 7^7* (i = 1, 2). 

In two dimensional space, the orbital angular momentum has only one component: 

L = e'^xY, (8) 



where e*-' is antisymmetric in i, j and e^^ = 1, and summation over the repeated Latin 
indices i, j is imphed. It is easy to show that [L,H] — ie^^a'^p' ^ 0, thus L is not a 
constant of motion even when = 0. Let 

S = -fH- (9) 

It is easy to show that [5", li\ = —ie^^a^p', then [L + S,H] = 0, and the quantity 

J^L + S (10) 

is a constant of motion. It is natural to regard J as the total angular momentum and S 
as the spin angular momentum in two dimensions. 

II. A. Solutions in the external central field 

To solve Eq.(6) a representation of the Dirac matrices is necessary. This can be reahzed 
by the Pauli matrices: 

70 = (7^ 7^ = i(7\ -f^^ia^. (11) 
In this representation S — Let 



\ G'(r)e^(^+V2)7^ 



It is easy to show that ipj is an eigenf unction of J: 



, ±1/2, ±3/2,.... 

(12) 

(13) 



So that ipj has angular momentum j. (There is no danger of confusing the angular 
momentum with the Latin index j used above, as the latter will not appear henceforth.) 
The radial wave functions F and G satisfy the system of equations 



F' - '—F +(E + ii- V)G = 0, 
r 



(14a) 



G' + —G -{E-ii- V)F = 0, 
r 



(14b) 



where j± = j ±1/2, and primes denote differentiation with respect to argument. We will 
have occasions to use the system of equations for vi and V2, so we also write down it here: 



v[--vi + {E + ii- V)v2 = 0, 



v!^ + ^-V2-{E- 1^- V)vi = 0. 



Prom Eq.(14a), we have 



G 



F' - {j_/r)F 



Substituting this into Eq.(14b) we get an equation for F alone: 



(15a) 

(15b) 

(16) 



F"+ - + 



V 



r E + n-V^ 



F' + 



j-V 



r(E + ii- V) 



F = 0. (17) 



With a given potential V{r) that is regular everywhere except possibly at r = and 
appropriate boundary conditions one can in principle solve Eq.(17) for F, and get G from 
Eq.(16). Equivalent to Eqs.(16) and (17), we may have 



G' + {u/r)G 
E-fi-V ' 



G' + 



E'-i,'-^-2EV + V' + - 

r{E — fj, — V) 



(18) 



G^O. (19) 



If in some region or at some point E + /j, — V — 0, we can not use Eqs.(16) and (17). But 
then E — ii — V^Om that region or in the neighbourhood of that point, and we can use 
Eqs.(18) and (19). Indeed, the fundamental equation is Eq. (14) or Eq. (15). They are 
regular everywhere except at r = 0. In principle one can solve them by direct integration 
without the help of Eqs. (16-19). Eqs. (16-17) or Eqs. (18-19) are to be employed when 
convenient. The possible singularities in these equations due to the vanishing denominator 
E + fjL — V or E — ijl — V should not cause any trouble in principle. Nevertheless, attention 



6 



should be paid to these possible singularities when we have to use Eq. (17) or Eq. (19) 
in the whole range of r. 

For free particles, V" = 0, then Eq.(17) becomes 

F" + ^F' + (^E^ - - F = 0. (20) 

In order to get well behaved solutions one should have E^ — fj,'^ > 0. Thus we have positive- 
energy solutions with E > /i and negative-energy solutions with E < —/i. Let us define 



k = \/E'^ — n'^ > 0, and denote positive- (negative-)energy solutions by the subscript k 
{—k), thus we have, say. 



E±k = ±Ek = ±^Jk^ + l^^. (21) 
It is not difficult to find the following solutions for free particles: 



.(0) /„N / (Ek ± l^)k 



nW^J^^^^Ukr), (22a) 



G%ir) = ±eij) f-^^^^Jm^Mkr), (22b) 

where m — Jm{kr) is the Bessel function, e(j) = 1(— 1) when j > 0(< 0), and the 
normalization factors are chosen such that the orthonormal relation takes the form 

/ drV'i'^,,(r)V'i%(r) = 6{k - k%,>, (23) 

Jdri;%{r)i;%{r) = 0. (23') 

The completeness of these solutions is ensured by the following relation which can be 
verified straightforwardly: 

E r dk [C(r)V'i?'(r') + ^%{r)^^-l]{r')] = S{v - v'). (24) 
When r — > oo, the radial wave functions have the asymptotic form: 



This can be obtained by using Eq.(22) and the asymptotic formula for the Bessel function: 

^"^^^™/^'°K^"T"l)- ^^^^ 

Now we consider particles moving in the external central potential V{r). We assume 
that V{r) — > more rapidly than when r — > oo, and is less singular than when 
r ^ 0. Then, for very large r, Eq.(17) takes the same form as Eq.(20). It is easy to see 
that E"^ > jji'^ gives scattering solutions while E"^ < ji"^ gives bound state solutions (bound 
states with E = ±/x are also possible, see Sec.V). Scattering states will be denoted as 
above, while bound states will be denoted by a subscript k which takes discrete values. 
The orthonormal relations are given by 

|dr4,,^.,(r)V^±,,(r) = 5{k - k%j,, (27) 

/drV^l,^.,(r)V^«,(r) = 5,,,5,,,, (27') 

and vanishing ones similar to Eq.(23'). The completeness relation is similar to Eq.(24) 
but has an additional term on the left-hand side (Ihs) : 

Y: r dk + ^-kAr)^Uj{r')] + E V'«.(r)V'l,(r') = S{r - r'). (28) 

As pointed out above, Eq.(17) takes the form of Eq.(20) at large r, so the solution 

Fkj{r) is giver by a linear combination of Jmikr) and Nm{kr), the Neumann function, at 
large r. Using Eq.(26) and 



and Eq.(16), the asymptotic forms for the radial wave functions can be shown to be 

(30a) 



'Ek±[x 



rriTT TT 



v±kji{r) \l cos kr- — --+ r]j{±Ek) 



2 4 



v±kj2[r) ±\l sm 



(30b) 



when r — > oo, where m — as before, rjj{±Ek) are the phase shifts. They depend on 
j rather than and also depend on the sign (not only the magnitude) of the energy, 
as Eq.(17) does. Compared with Eq.(25), the asymptotic forms in the external certral 
field are distorted by the phase shifts. But it should be remarked that the normalization 
factors in Eq.(30) are the same as in Eq.(25). 

II. B. Partial-wave analysis of scattering by the central field 

It is well known that positive-(negative-)energy solutions correspond to particles (an- 
tiparticles) after second quantization. In this subsection we discuss the scattering of 
positive-energy solutions by the central field V{r) described above. The scattering of 
negative-energy solutions can be formally discussed in a similar way. 
The probabihty current density associated wiht Eq.(2) is given by 



(31) 



The incident wave may be chosen as 



V 2E, 



A. 



6^^=^, A; > 0, 



\ \l 2Ek J 

which is a solution of Eq.(6) with positive energy E^ when x 
probability current density is 

jm ^xk/ Ef; GxV, 



(32) 



-oo. The incident 



(33) 



where is the unit vector in the x direction, the incident direction in this case. The 
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(34) 



scattered wave should have the asymptotic form when r — > oo: 

^ " V f2{0) J 

where the factor y/i — e^'^^^ is introduced for latter convenience. Then the r and 9 
components of jsc at large r can be shown to be 



(35) 



However, the expression (34) does not solve the Dirac equation at large r if fi{9) and 
f2{d) are independent of each other. In order to satisfy Eq.(6) at large r, one should have 



f2iO) 



ik 



(36) 



Setting 



it is easy to show that 



jsc^erv\f{9)\yr 



(37) 



(38) 



when r — > oo, where is the unit vector in the radial direction. Thus the differential 
cross section (in two dimensions the cross section may be more appropriately called cross 
width) is given by 

a{e) = \f{e)\'. (39) 

The outgoing wave comprises Eqs.(32) and (34). With the relation (36) and the definition 
(37), it takes the following asymptotic form when r — > oo: 



2 Eh 



'le^'=V(^) 



2Ek 



(40) 
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On the other hand, the solution of Eq. (6) with definite energy Ek > has the form 

^|;{r,e) = J2aj^kj{r,e), (41) 
j 

where the summation is made over all j = ±1/2, ±3/2, The asymptotic form of 

Eq.(41) can be obtained from Eq.(30). By using the formula 

^ikx^ J2 il"lj|„|(A;r)e^"^ (42) 

n=— oo 

and Eq.(26), one can compare Eq.(40) with the asymptotic form of Eq.(41). They must 
coincide with each other for appropriately chosen a/s. In this way one finds all Uj in 

terms oir]j{Ek) and 

/(^) = E V^^'"'^"''^^ nAEk)e'^-'- (43) 
The total cross section at turns out to be 



at = J^ dea{0)^-J2sm'VjiEk). (44) 



^27r _ 4 

One easily realizes that all information of the scattering process is contained in the phase 
shifts. The latter are determined by solving the system of equations (15) with the bound- 
ary conditions (30). The purpose of the Levinson theorem is to establish a relation between 
scattering states and bound states, specifically, to establish a relation between the phase 
shifts and the total number of bound states in a given angular momentum channel j. 

III. Phase shifts near threshold 

In this section we discuss the behavior of the phase shifts r)j{±Ei;) near A; = 0. This will 
be employed in the next section. For exact analysis let us cut off the potential. That is, 
we consider potentials that satisfy V{r) = when r > a > 0. Such potentials will be 
denoted by Va{r) in the following. In the region r > a, then, Eq.(17) reduces to the form 
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of Eq.(20), and the solution may take the form 



Flkj{r) = ^^^||^[cosr;,(±E,) J^(fcr) - sinr^,(±^,)7V„(A;r)], (45) 

where the superscript ">" indicates r > a. Using Eqs.(16), (26), and (29), one can obtain 
the expected asymptotic forms (30). In the region r < a, Eq.(17) cannot be simphfied. 
Let us consider the behavior of the solution near r = 0. We have assumed that V{r) is 
less singular than when r — > 0. Thus, when r — > 0, V{r) is regular or behaves like 
Uo/r^ where Uq is a constant and < 5 < 1. Accordingly, Eq.(17) becomes to leading 
terms 

1 ?2 

F" + -F' - -^F = (46a) 

in the first case or 

F" + i±^F' - '-^'-/^^F = (46b) 

in the second case. Therefore, the regular solution of Eq.(17) may have the following 
power dependence on r when r — > 0: 

/f(r,A;)^r"^ (47a) 

in the first case or 



(47b) 



j < 



in the second case, where we have denoted F±kj{r) with these boundary conditions by 
f^{r, k), and m — as before. The solution of Eq.(17) in the region r < a is 

F^kjir)=Af{k)f^{r,k), (48) 

where the superscript "<" indicates r < a. In general the coefficient Af depends on k 
such that the two parts of F±kj{r) can be appropriately connected at r = a. Obviously, 
F±kj{^) and G±fcj(r) should be continuous at r = a, so that the probability density and 
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the probability current density are continuous at r = a. For simplicity, we assume that 
Va{r) is continuous at r = a, which means Ki(r) — > when r — > a~ and V^(a) = 0. Then 
from Eq.(14) we see that is also continuous at r = a. Therefore, -f±fcj(^)/-^±fci(^) 

is continuous at r = a. This leads to 

where ^ — ka and 

m = f^. (50) 

where the prime indicates differentiation with respect to r. The above result shows that 
the behavior of rjj{±Ef;) is determined by that of Pf{^) and ultimately by that of (a, k). 
The general dependence of ff{r, k) on k may be very complicated since Eq.(17) depends 
on /c in a rather complicated way. This is quite different from the nonrelativistic case 
where the function fm{i^,k), the counterpart of f^{r,k), is an integral function of k[16]. 
Fortunately, only the property of (r, k) near A; = is necessary for our purpose. 

We consider the limit E ^ /j, oi Eq.(17). In this limit it takes the following form to 
the first order in k'^: 

F" + P{r,k)F' + Q{r,k)F = 0, (51) 



where 



Q{r, A;) = - 2pLV + V' - 



1 



fi 2^r(2^ - V) 



k\ (51'b) 



r2 ^ r(2/i - V) 

We denote the solution of this equation that satisfies the boundary condition (47) by 
fj'ij-, k). Note that both P(r, k) and Q{r, k) are integral functions of /c, and the boundary 
condition (47) is independent of k. Then a theorem of Poincare tells us that fj~{r,k) is 
also an integral function of k. On the other hand, Eq.(17) coincides with Eq.(51) in the 
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limit E —>■ II. Therefore, f^{r,k) must coincide with fj'{r,k) in the hmit A; — > since 
they satisfy the same boundary condition. Hence we conclude that //(r, k) is an analytic 
function of k in the neighbourhood of A; = 0. Moreover, fj^{r, k) is an even functions of 
A;, since Eq.(17) is invariant under the change k ^ —k, and the boundary condition is 
independent of k. The above conclusion holds regardless of whether V{r) = K(r) or not. 
We can now proceed as in the nonrelativistic case [16] and arrive at the result 

te.nrjj{Ek)^b+e'^ or ^ (^^0), (52a) 

where 6^ 7^ is a contant and is a natural number. For a strong repulsive potential 
Eq. (51) may become singular at the points where V — 2/1. Then the above analysis is 
not reasonable. In this case, however, we may consider Gkj{r) instead of Fkj{r) and study 
the limit E ^ ^ oi Eq. (19). The same result (52a) can be attained in a similar way. A 
similar analysis leads to 

tanr/,(-Efc)^67C'"7 or ^ (^^0), (52b) 

where 6j 7^ is a contant and is a natural number. 

In the above analysis we have assumed that Va{r) is continuous at r = a. This 
is, however, not mandatory. Using the condition that both F±kj{r) and G±kj{r) are 
continuous at r = a, one can arrive at the same result (52) even when V{r) is not 
continuous at r = a. 

IV. The Levinson theorem 

With the above preparations, we now proceed to establish the Levinson theorem by the 
Green function method. We define the retarded Green function G(r, r', E) of the Dirac 
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equation in an external potential V{r) by 

G(ry.E) = Y.^f^^. (53) 

where {V'rlr)} is a complete set of orthonormal solutions to Eq.(6) where V may be a 
noncentral potential, Er is the energy eigenvalue associated with the solution '0T(r), and 
6 = 0"*". G{r,r',E) satisfies the equation 

{E-H + ie)G{r,r',E)^5{r-r'). (54) 

For a free particle, V — 0. We denote the Hamiltonian by Hq and the solutions to the 
free Dirac equation by ip!^^ (r) . The retarded Green function in this case is defined as 

where E!^^^ is the energy corresponding to the solution ■0|°^(r). G(°)(r, r', E) satisfies 

(E-Ho + ie)G^'\r, r', E) = 5(r - r'). (56) 

We have the integral equation for ^(r, r', E): 

G{r,r',E)-G^''\r,r',E) = J dr" G'^°\r,r" , E)V{r")G{r" ,r' , E). (57) 

In a central field V^(r) = V{r) [not necessarily Kj(r)], we have discussed the solutions 
■0r(r) of the Dirac equation (6) in Sec.II. They contain two different classes: ip±kj{r,0) 
with continuous energy spectrum E±f; and ipf^j{r,9) with discrete energy spectrum E^j. 
However, we may discretize the continuous part of the spectrum by requiring the radial 
wave function E±kj{T) or v±kji{r) to vanish at a sufficiently large radius [i? ^ a for 
T4('^)]- In this case we will denote all the solutions by t/j^jir^O) and the corresponding 
energies by E^j. We also apply the same prescription to the free Dirac particle to turn 
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the continuous energies E^^l into discrete ones 

E^S- The waves functions will be denoted 

by (f, 0) in the free case. We have then 

/ G,n(r, r', E) G,i2(r, r', E)e-'<^' \ 

G(r,r',E)=;^^ , (58) 

\Gi2i{ry,E)e'' Gj22(r,r',E)e'('-''^ J 

where 

^^^^ ^ ^ V^{E-E,j + zey ^ ' 

where p, g = 1,2 are used to denote spinor indices. For a free particle we have similar 
relations. The integral equation for Gjpq{r, r', E) may be derived from Eq.(57). The result 
turns out to be 

GUr.r'.E) - G^S,(r,r',E) = / dr" r"^G'Ss{r,r",E)V(r")Gjs,{r"y,E). (59) 
Since the wave function should be finite at r = 0, we have the boundary conditions 

i;«,,(0) = 0, t;Si(0)=0. (60) 

For bound states with \Enj\ < /x, we have Vnjp{oo) = 0. For scattering states with 
{E^jl > II we may set v^jpir) — when r > R. However, it should be remarked that in 
general ^^^^(-R) 7^ when we require v^ji{R) = 0. In any case, we have 

v^jpioo) = 0, v'^SpM = 0- (60') 

Now the orthonormal relation takes the form (27') for both bound and scattering states, 
from which we have 

{v^>j,v^j) ^ S^^', Vj, (61) 

where 

poo 

{v^,j,v^j)= dr^vl,jp{r)v^jp{r). (62) 

JO rt 
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As we have set v^jpir) = when r > i? for scattering states, the upper bound of integration 
in Eq.(62) need not be replaced by it! for these states. For the free particle we have 

= Vj. (63) 

Using Eqs.(58') and (61), it is easy to show that 

fdrrY: G,pp{r, r,E)^Y. ^ ^ ^ , ■ ^^^^ 



Employing the mathematical formula 



^-^P--m5{x), (65) 



X -\-ie X 

and taking the imaginary part of the above equation, we have 

Im [drrY^ Gjpp{r, r, E) = -ttJ^ S{E - E,^). (66) 

Integrating this equation over E from — /x~ to /x" yields 

Im / dE I drr^^Gjpp{r,r, E) = —njir, (67) 

where nj is the number of bound states in the angular momentum channel j with \Ef^j\ < 
II, li~ — II — 0, and —iJ,~ — —II + 0. The existence of bound states with critical energies 
E^j = ±11 does not alter this result. Similarly, we can show that 

Im r dE [ drrY,GfJp{r,r,E) = 0. (68) 
J-n- J p 

Combining Eqs.(67) and (68) we obtain 

Im r dE fdrr Y.[Gjw{r, r, E) - G^^p{r, r, E)] = -njii. (69) 
J-n J p 

On the other hand, substituting Eq.(58') and a similar one for the free particle into the 
right-hand side (rhs) of Eq.(59) we get 

/ drrY.lG„[r.r.E) - GZ(r,r,E)] = £ ''g;^"^'';"";^'' . (70) 

•f p [E - + te) [E - E^,j + le) 
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where the definition of {vi^ij^v).-) is similar to Eq.(62) and 

{v^'hVv^'i) = r^'' EC(^)^(^)^'^Op(r)- (71) 

JO p 

Using the system of equations (15) and a similar one for the free case, and employing the 
boundary conditions (60) and (60') , one can show that 

(^S\ Vv^,,) = {E^,, - 4?)(^g\ v^,,). (72) 

Substituting this result into Eq.(70) and taking the imaginary part, we get 

P kk' 

(73) 

Integrating this equation over E from — oo to +oo it turns out that 

Im dE jdrr E[G,pp(r, r, E) - CfJ^ir, r, E)] = 0. (74) 

J —oo J p 

As in the nonrelativistic case, this means that the total number of states in a specific 
angular momentum channel is not altered by an external field, except that some scattering 
states are "pulled down" into the bound-state region. Here the external field may be either 
attractive or repulsive. This is different from the nonrelativistic case, where bound states 
exist only in attractive fields. Another new character of Dirac particles is that when 
the potential, say, an attractive potential, becomes strong enough, some bound states 
may even be pulled down into the region of negative-energy scattering states. This new 
character, however, does not alter the above conclusion that the total number of states 
in a specific angular momentum channel remains unchanged, since the conclusion is a 
consequence of the Dirac equation (15) and the completeness of the whole set of states. 
If the integration over E in Eq.(74) is performed from — to /i", we have 

Im r dE I drrY}G,^{r,r,E)-Gfl{r,r,E)] = -vr^'E IK'i'^i?)!'' (75) 
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where we have taken into account the fact that {Ej^h > /x, and the prime in YI'k' indicates 



that the summation over k' is performed only for bound states with \Ei^ij\ < /i. Using the 
completeness relation (28), which now takes the form 

EV'«.WUr')=^(r-r'), (28') 
where the summation over k is performed for all states, we can show that 

,(0) /'^^.,(o)*,'^'^ 



E 



This in turn leads to 



(76) 



(77) 



Substituting Eq.(77) into Eq.(75) we recover Eq.(69). Combining Eqs.(69) and (74) we 
obtain 





r-fj, r+oo 




Im 


y +7- 


iEj 




J —oo J n 





dE I drr YliGjppir, r, E) - CfJ^ir, r, E)] = njTr. 



(78) 



We have thereupon finished the first step in our establishment of the Levinson theorem. 

The next step is to calculate the Ihs of Eq.(78) in another way. In the above treatment 
we have discretized the continuous spectrum of E^^^ and the continuous part of Ei^j. In 
the following we will directly deal with these continuous spectra, and use the notations 
of Sec.II. Then the retarded Green function G{v,v',E) is given by Eq.(58) but where 
Gjpq{r, r', E) is given by 



+ E 



dk 



+ 



v-kjp{r)v*_kjq{r') 
_V^{E - Ekj + ie) ' V^{E - E^kj + ie) 
VKjp{r)v*jq{r') 



(79) 



where the integration is performed over scattering states while the summation over bound 
states, and E±kj — E±k [cf. Eq.(21)] is independent of j. For Gfjq{r,r', E) we have a 
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similar expression but without the last term. Using the formula (65) it is easy to show 
that 

Im / drr^Gjpp{r,r,E) 
■I p 

roo 

= -TT / dk [5{E - Ekj){vkj, Vkj) + S{E- E_kj){v-kj, v^kj)] -7iJ2S{E- E^j). (80) 

Here the inner products {vkj, Vkj) etc. are defined in the same way as Eq.(62). Integrating 
this equation over E as follows we have 



Im 



J — oo J n 



dE / drr^Gjpp{r,r,E) 
■' p 

roo 

-n / dk [{vkj, Vkj) + {v_kj,v-kj)] - 7r(5^j + 5_^j), (81) 

t/ 



where S/^j = 1 (5-^j = 1) if there exists a bound state in the angular momentum channel 
j with critical energy E^j — n {E^j — — //), otherwise — {S-^j — 0). (See Sec.V for 
further discussions.) Similar to Eq.(81) we have for the free case 

+ / dE drrY.G'S^{r,r,E) = -tt / dkHv'^^ ,v'^^) + {v%,v%)]. (82) 

i_ -oo J IJ, i p ^ 

Combining Eqs.(81) and (82) we obtain 

rY.[Gj,,ir,r,E)-GfJ^{r,r,E)] 

p 

lo '"'"^'^ ~ 4?)] - ^ dk [{v^kj, v-kj) - {v%, v%)] 

-7i{d^,+S_^,). (83) 

From the orthonormal relations (23) and (27) one can show that 

{v±k'j,v±kj)^S{k-k'), {v%i,v%) = 5{k-k'), Vj. (84) 

This implies that both integrands in the rhs of Eq.(83) are 5(0) — 5(0). There is, however, 
a subtle difference between these two 5(0) 's, and it is this difference that leads to the 
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r—^i r+oo 




Im 


y +7- 


iEj 




J —oo J ^ 





Levinson theorem. To get rid of the difficulty of infiniteness, we define 

{vij,Vkj)ro= dr22'"ijp(^Hjp{r), (85) 

JO p 

and obtain {vkj,Vkj) in the hmit / — >• k and ro — > cxd. Using two systems of equations 
satisfied by Vkjp and vijp [cf. Eq.(15)], and the boundary conditions [cf. Eq.(60)] 

Vkj,{0)^0, 4>)=0, (86) 

it can be shown that 

{Eij - Ekj){vij, Vkj)ro = v!j2{ro)vkji{ro) - vlji{ro)vkj2{ro) ■ (87) 
Since Tq is large, we can use Eq.(30) to evaluate the rhs and in the limit Z — > A; we get 

{vkj, v,,)r, = ^ + l^^^ - i-r^^ co«[2^^o + 2r,,{E,)]. (88) 
Similarly, we have 

(4°',4°'k = ^-(-r^cos2*., (89) 

Obviously, the infiniteness lies in the first term To/tt when we take the limit Tq ^ oo. 
This disappears when we subtract Eq.(89) from Eq.(88). Using the well-known formulas 

sin2A;ro , , 

lim = 6(k), 90 

and g{k)5{k) = g{0)d{k) for any continuous function g{k), we obtain 

= 1^!?^ + tlls{k) sin277,(/x) + (-)™^ cos2A;rosin2r/,(i?,). (91) 

TT CLK Z TTKrjk 

In a similar way, we can show that 

/ ^ _ r (°) ^ 

[V-kj,V-kj)ro {V_kj,V_i^j)ro 

= _ t}ls{k)sm2rjj{-pi) - (-)m^cos2krosm'rjj{-E,). (92) 
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So far in this section V{r) need not be Va{r). In the following we set V{r) = V^(r). Then 
Eq.(52) holds. Therefore sin 277j(//) = sin 277j(— //) = 0. Integrating Eq.(91) over k (from 
to +00) and taking the limit ro ^ 00 we have 

/ dk[{vkj,Vkj) - {v^kj,Vkj)] 

J U 

= ^fe(+oo) -r^,(/.)] + {-r^^ hrn^f dk '^'^^[^'^ cos2A;ro. (93) 

The last term in this equation can be decomposed into two integrals, the first from to 
£ = O'*', while the second from e to +00. The second integral vanishes in the limit ro — > 00 
since the factor cos 2/cro oscillates very rapidly. For the first integral, wc have 

sin^ iljiEk) _1 r sill- iij(Ek) _ I f^%. sin^ I ij(Ek) 



rdfe"">^^^^os2fcro=i Tdfe'^'f ^^-^ = i rde. , 

Jo kEk jJ, Jo k jj, Jo ^ 

as A; is very small. For the same reason we can use Eq. (52a). In the first case, iwir]j{E}^) — 
we have s=ni^r]j{Ek) = bf^^t, so that 

I .-^^si^.^^ ^ 6f ^ ill . ^ ^_ 

II Jo ^ II Jo 4:npj 

In the second case, t&nrij{Ek) = 7r/21n^, we have sin^ rij{Ek) = 7r^/41n^^, so that 

sin^^^,™ 1 1 ^ 

H Jo ^ An Jo ?ln ? Aninea 

where the final integral is calculated by setting In^ — w. Thus the first integral vanishes 
as well. Therefore, we have 

dk[{v,„Vkj) - (4°\4°))] = kvA+oo)-vAl^)]- (94) 



On the basis of Eqs.(92) and (52b) we can show that 

^ dk [{v-kj, v-kj) - {v%, v%)] = ^[r]j{-Qo) - r^j(-//)]. (95) 
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Substituting Eqs.(94) and (95) into Eq.(83) we obtain 

dE I drr YliGjppir, r, E) - CfJ^ir, r, E)] 





r — li r+oo 




Im 




cIeJ 




J —oo J 1^ 





p 

= hifJ') - Vj{+<^)] + hi-fJ') - - n{S^j + 5^^j). (96) 

Combining this result with Eq.(78) we arrive at 

hilJ') - Vj{+oo)] + [Vji-fJ') - ?7j(-oo)] = rijn, (97) 

where 

Uj = nj + Sf,j + S-^j (98) 

is the total number of bound states, including the possible ones with critical energies ±/x, 
in the angular momentum channel j. Eq.(97) is the Levinson theorem for Dirac particles 
in an external central field in two dimensions. In the next section we will discuss some 
relevant points of the theorem. 

V. Discussions 

1. On bound states with critical energies. We consider the cut-off potential Va{r). In the 
region r > a, \4(r) — 0. Equations (16) and (17) are applicable for E — /i, while for 
E — —fj, Eqs.(18) and (19) should be employed. 

In the region r > a and for E — /i, Eq.(17) becomes 

F" + If' - = 0, (99) 

where m — | j_ | as before. The well behaved solution is 

F>(r) = r-"^. (100) 

The corresponding solution G^j{r) obtained from Eq.(16) decreases more rapidly than 
F^j{r) when r — > oo. Thus the normalizability of the solution is determined by the 
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behavior of Fj^j{r) at large r. Obviously, the solution can be normalized only when 
m > 1. It is well behaved when m — 0,1, but cannot be normalized. In other words, the 
solution with energy E = n is a. bound state only when j > 3/2 or j < —1/2. The cases 
j — ±1/2, 3/2 do not correspond to bound states. Of course, a bound state with E — ji 
actually exists only when the potential Va{r) has a specific form such that the solution is 
continuous at r = a. For a square- well potential with depth Vq, this leads to 

2m + Jmikoo) (101a) 

for j > 3/2 or 

Jm-i{ha) = (101b) 



for j < —1/2, where ko = \JVq + 2/iVo. For a given j, this is satisfied only for some 
specified depth Vq or radius a. 

For E = —fi, we consider Eq.(19). In the region r > a it reduces to 

1 m'^ 

G" + =-G' - -^G = 0, (102) 

where m' — \ The well behaved solution is 

Gl,,{r)^T--'. (103) 

One can get F^^j{r) from Eq.(18). It decreases more rapidly than G^^j{r) when r — > oo. 
Therefore the solution is normalizable when m' > 1. In other words, the solution is a 
bound state when j > 1/2 or j < —3/2. The cases j — ±1/2, —3/2 do not correspond 
to bound states. For the square-well potential, a bound state with E — —/i really exists 
only when Vq > 2/1 and satisfies 

koaJm'+i{koa) = (^m' - Jm'iha) (104a) 
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for j < —3/2 or 

Jm'-lM = (104b) 

for J > 1/2, where kg — \JVq — 2/xVo. Given j, say j — —5/2, there exist infinitely many 
solutions of Vo to Eq. (101b). They are functions of a. In general most of them cannot 
satisfy Eq. (104a) at the same time as the two equations are independent, thus the two 
critical energy bound states do not appear simultaneously. However, the solutions of Eq. 
(104a) are also functions of a. By varying the parameter a, it may be possible to match 
some specific solution of Eq. (101b) with some specific one of Eq. (104a). When this 
really happens, the two critical energy bound states can appear simultaneously. Anyway, 
it cannot be asserted that the two critical energy bound states never appear simultaneously 
for any potential. 

The form of the Levinson theorem (97) is not modified by the existence of the critical 
energy states, regardless of whether they are bound states or not. The existence of critical 
energy bound states just changes Uj from nj to + 1 (when there is one with E — fx 
or E — —/i) or nj + 2 (when the two appear simultaneously), and does not alter the 
form of Eq. (97). In three dimensions, the theorem involves additional terms that vanish 
except when there exist half bound states[10]. In two dimensions there is no such term, 
which is clear from the result (97). Thus the theorem in two dimensions is not affected 
by the existence of half bound states. The reason for the difference between two- and 
three-dimensional cases is similar to that in the nonrelativistic theory. This has been 
discussed in detail in Ref. [16]. 

2. About r7j(±oo). We write down two systems of equations for the radial wave 
functions v^j Jr) and v^jp{r) in two external fields U{r) and U{r), respectively. Using the 
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boundary condition (86) and the asymptotic form (30), it is not difficult to show that 

sm[fjf{E,) - rj^iE,)] = -^(^S, {U - U)v^^). (105) 

Here the superscript U is used to distinguish the wave functions and phase shifts from 
those in the external field V{r). Now we set U{r) = XV (r) [V{r) not necessarily be 
Va{r)], U{r) = (A + AX)V{r), and denote vgp{r) = Vkjp{r, A), r}f{Ek) = r]j{Ek, A). When 
AA — > 0, v^jp{r) = Vkjp{r, A + AA) can be replaced by Vkjp{r, A), and Eq.(105) becomes 

sm[A7jj{E,, A)] = -^AA(^;fe,(r, A), V{r)vkj{r, A)), (106) 

where Ar]j{Ek, A) = r)j{Ek, A + AA) -r)j{Ek, A). Obviously, Vkjp{r, 0) = Vkjp{r, 1) = 

Vkjp{r), r)j{Ek,l) — r]j{Ek). It is natural to define r)j{Ek,0) = in the absence of an 
external field. It is also natural to require that for any finite k, rij{Ek,X) be continuous 
functions of A when A varies continuously from to 1, as all quantities in the equation 
and boundary conditions are continuous in A. With this requirement and the above 
definition r]j{EkiQ) — one can determine r]j{Ek). It should be remarked, however, that 
r]j{ii, A) = limfc^o^j(-Sfc) ^) is not continuous in A. Otherwise the Levinson theorem would 
be impossible. To see this, let us have a look at the original Levinson theorem in an 
attractive potential XV{r) in three dimensions: 

(5i(0,A)-5i(oo,A) =n,(A)7r, Z = 0,1,2,..., (107) 

where we have not included the modified I — case. If (5^(0, A) is a continuous function of 
A, then the Ihs of Eq.(107) is continuous in A. On the other hand, n;(A) is a nonnegative 
integer. When A varies continuously from to 1, n;(A) varies from to ni by discontin- 
uous jumps. The rhs of Eq.(107) is obviously not a continuous function of A. This is a 
contradiction. So that (5;(0, A) cannot be continuous in A. The case for r]j{^, A) is similar. 
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Since r]j{Ek, A) is a continuous function of A, the sm[Ar}j{Ek, A)] on the Ihs of Eq.(106) 
can be safely replaced by ArjjlEk, A), and we have 

dr]j{Ek,X) _ nEk 



dX k 



'ivkj{rA),V{r)vki{r,X)). (108) 



When Ek is infinitely large, since V{r) is not very singular at r = and is regular 
elsewhere, we may ignore XV {r) in the system of radial equations and approximately 
replace Vkj{r, A) by vf^{r) in the above equation. Then Eq.(108) can be easily integrated 
over A from to 1 and results in 

= -^(4°>.V«g*). (109) 

Substituting the exact solution (22) into Eq.(109) and then replacing the Bessel functions 
by their asymptotic forms as k is large, we arrive at 

/•oo 

ry.(+oo) = -/ drV{r). (110a) 
In a similar way we can show that 

/■oo 

r]j{-oo) = / drV{r). (110b) 

J 

As we have assumed that V{r) is less singular than r^^ when r — > and decreases more 
rapidly than when r — > oo, the above integrals converge. These results have the same 
form as those in three dimensions [9, 17, 18]. Of course, Eq.(llO) holds in the special case 
V{r) = Va{r). As a consequence of the above results, we have 

r^,-(+oo)+7^,-(-oo) =0. (Ill) 

This reduces the Levinson theorem (97) to the form of Eq.(l). It means that the sum 
of the phase shifts at the two thresholds serves as a counter for the bound states in a 
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specific angular momentum channel. This is similar to the case in three dimensions, but 
is somewhat simpler. 

The Levinson theorem for Dirac particles could not be separated into two parts, each 
of which hkes that for Schrodinger particles [6,9,10,19]. Essentially this is because that 
positive-energy solutions or negative-energy solutions alone do not form a complete set. 
An evidence can be seen as follows. In an attractive potential, say, Eq. (108) shows that 
r)j{Ek,X) is positive and increases with A. Similarly, it can be shown that rij{—Ek,\) 
is negative and decreases when A increases. Thus r]j{—^) may be negative when the 
attractive potential becomes strong enough. In three dimensions this has been verified 
by numerical calculations [10, 18]. As r]j{—^) may be negative, it cannot always equal a 
nonnegative integer (in unit of tt) and thus cannot serve as a counter. 

3. Extension to more general potentials. Throughout this paper we have assumed that 
V{r) is less singular than r~^ when r — > and decreases more rapidly than when 
r — > oo. In the development of the Levinson theorem, we further cut off V{r) to the 
special case Va{r) for the sake of exact analysis. However, the radius a beyond which 
Va{r) vanishes is not specified in our discussion. Though both sides of Eq.(l) depend 
on the particular form of Va(r) and thus depend on a, the equality between them does 
not. Thus we expect that Eq.(l) remains valid when T4(r) is varied continuously to the 
limit a — > oo provided that the asymptotic form (30) holds and Uj remains finite in the 
process. This requires that V{r) decreases rapidly enough when r — > oo. It seems that 
the Levinson theorem holds at least for short-range potentials that decrease more rapidly 
than when r — > oo. 
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